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violating interactions is studied in the presence of the gauge mediated super- 



^^ ■ symmetry breaking. Neutrino masses arise in this model at tree level through 

the RG-induced vacuum expectation values of the sneutrinos and also through 
direct contribution at 1-loop. The relative importance of these contributions 
is determined by the values of the strong and weak coupling constants. In 



^ ■ case of purely A' couplings, the tree contribution dominates over the 1-loop 

diagram. In this case, one simultaneously obtains atmospheric neutrino oscil- 
lations and quasi- vacuum oscillations of the solar neutrinos if all the A' cou- 
plings are assumed to be of similar magnitudes. If R parity violation arises 
from the trilinear A couplings, then the loop induced contribution dominates 
over the tree level. One cannot simultaneously explain the solar and atmo- 
spheric deficit in this case if all the A couplings are of similar magnitude. This 
however becomes possible with hierarchical A and we give a specific example 
of this. 



I. INTRODUCTION 

Variety of observations of the solar and atmospheric neutrinos have given important 
information on the possible structures of neutrino masses and mixings. Based on these 
observations, one can infer that oscillations among three active neutrinos are likely to be 
responsible for the observed features of the data. These oscillations are characterized by 
two widely separated mass scales. At least one mixing angle involved in oscillations of at- 
mospheric neutrinos is large [ij. The detailed data on the day-night asymmetry and recoil 
energy spectrum of the solar neutrinos seem to favour the presence of one more large mix- 
ing angle 0. Thus the neutrino mass spectrum seems to be characterized by hierarchical 
(mass) 2 differences and by two large and one small mixing angles, with the small mixing 
angle demanded by the CHOOZ experiment ||. Many mechanisms have been advanced to 
understand these features of the neutrino spectrum || . One of these is provided by super- 
symmetric theory which contains several features to make it attractive for the description of 
the neutrino spectrum. The lepton number violation needed to understand neutrino masses 
is in-built in this theory through the presence of the R parity violating couplings |J . More- 
over, it is possible to understand the hierarchical neutrino masses and large mixing among 
them within this framework without fine tuning of parameters or without postulating ad 
hoc textures for the neutrino mass matrices ||[7|]. 

The supersymmetrized version of the standard model contains the following lepton num- 
ber violating couplings: 

Wf, = €iLiH 2 + \ijkLiLje c k + \ i j k L i Qjd'j s ., (1) 

where L, Q, H 2 represent the leptonic, quark and one of the Higgs doublets (up-type) respec- 
tively and e c , d c represent the leptonic and down quark singlets. Each of these couplings is a 
potential source for neutrino masses. There have been detailed studies of the effects of these 
couplings on neutrino masses under different assumptions @-ll|. We briefly recapitulate 



the relevant gross features of these studies and motivate additional work that we are going 
to present in this paper. 

The most studied effect is that of the three bilinear mass parameters q, particularly in the 
context of the supersymmetry breaking with universal boundary conditions at a high scale 
P-P| JiT||rT| . This formalism provides a nice way of understanding suppression in neutrino 
mass m v compared to the weak scale. The mixing among neutrinos is largely controlled in 
this approach by the ratios of the parameters e^. This can be approximately described by 
the following matrix for a large range in parameter space of the minimal supersymmetric 



standard model with radiative SU(2) x U(l) breaking |6|,|7|.pl| : 



U 



( \ 

C\ S±C 2 S1S2 

— Si C\C 2 C\S 2 

— S 2 C 2 



(2) 



where 



si 



ei 



4 + 4 



S-2 



4 + 4 



\ 4 + 4 + e. 



2 • 
3 



(3) 



The above matrix can nicely reproduce the small angle MSW solution together with the 
atmospheric neutrino anomaly. However, due to the specific structure, one cannot have two 
large mixing angles keeping at the same time \U e s\ < 0.1 as required from the negative results 
of CHOOZ. Thus purely bilinear R violation with universal boundary conditions cannot 
account for the observed features with two large mixing angles. This may be remedied 



by not insisting upon the universal boundary conditions [ 10 1 or by adding right handed 



neutrinos 12 . 



In contrast to the bilinear case, the presence of trilinear interactions can allow two 
large angles without conflicting with the CHOOZ result. There have been several studies 
to determine possible set of trilinear couplings which can reproduce the observed features 



of neutrino masses and mixing |L3]. It is not surprising that one could 'fit' the neutrino 
spectrum in these cases due to very large number of trilinear couplings. But it was realized 
|T4] , |T5|| that gross features of the neutrino spectrum can be understood without making 
specific assumptions on the trilinear couplings other than requiring them to be similar in 
magnitude. This makes the R violation with trilinear interaction 'predictive' in spite of the 
presence of very large number of couplings. 

In addition to the R violating parameters, the neutrino spectrum in these models also 
depends upon the nature of the supersymmetry breaking. This spectrum has been studied 
in the standard supergravity case with bilinear as well as trilinear couplings and in the case 
of gauge mediated supersymmetry breaking when R violation is only through the bilinear 
terms in eq. ([!]). The supersymmetry breaking generically introduces two different types of 
contributions to neutrino masses. The presence of terms linear in the sneutrino field in the 
scalar potential induces a vacuum expectation value (vev) for the former which mix neutrinos 
with neutralinos and lead to neutrino masses. In addition to this 'tree level' contribution, 
the trilinear terms in the superpotential also lead to neutrino masses at the 1-loop level. 



Both the types of contributions are present in models with purely bilinear or purely trilinear 
terms in the superpotential at a high scale. In case of the trilinear couplings, the running of 
the couplings to lower scale generates in the scalar potential couplings linear in the sneutrino 
vev and lead to a tree level contribution which is often neglected in the literature. 

The relative importance of the tree level and loop induced contributions to neutrino 
masses in case of trilinear interactions was studied in |15|,|16j in the context of the standard 
supergravity (mSUGRA) models with universal boundary conditions at a high scale. It was 
concluded that the tree level contribution dominates over the loop for large ranges in the 
parameters of the model. This results |T5[ in the following hierarchy in neutrino masses if 



all the trilinear A' couplings are assumed to be similar in magnitudes: 



TTTiA, Wlnnn TTh x 



U'z "'loop 



(4) 
m U3 m + m loop m b 

The parameters m loop and m® characterize the strength of the 1-loop and the RG-induced 

tree level contributions respectively and are determined by soft SUSY breaking parameters. 

m s and m& denote the strange quark and the bottom quark masses respectively. It is possible 

to obtain the vacuum or the MSW solution (large angle) to solar neutrino problem in this 



context by choosing the SUSY parameters in appropriate range [15] . 



An attractive alternative to the standard supergravity induced SUSY breaking is pro- 



vided by the gauge mediated SUSY breaking |21| . Neutrino mass spectrum has been studied 



in gauge mediated models with trilinear R- violation by Choi et al. in p2| . Their study has 
been confined to non-minimal models of this category. The minimal model in this category 
called the Minimal Messenger Model(MMM) E3| has only two free parameters and is more 



predictive than the standard SUGRA based models and the models studies in |22|]. The 
two free parameters of the model determine all the soft terms at the high scale ~ a few 
hundred TeV, where SUSY breaking occurs. Thus this model implies very constrained spec- 
trum for neutrino masses. This constrained spectrum has been shown to be inadequate 
for simultaneous solution of the solar and atmospheric neutrino anomalies in case of purely 
bilinear R violation. In this work, we wish to study neutrino masses in the minimal messen- 
ger model in the presence of purely trilinear R violation. This would mean that both the 
scale at which SUSY breaking occurs as well as the boundary conditions at the high scale 
would be sufficiently different from the mSUGRA scenario which has been studied in |L5| . 



We have studied the neutrino mass spectrum in the MMM for two separate cases, namely 
purely A' couplings and purely A couplings. Such a choice has been made for simplicity. The 
A' couplings with comparable magnitude are argued to describe neutrino spectrum well. In 
contrast, we find that if all the A couplings are of similar strength, then one cannot describe 



the neutrino spectrum well and one needs to postulate somewhat inverse hierarchy among 
them. We give a specific example with hierarchical A which reproduces the observed features 
of the neutrino spectrum. 

Within the Minimal Messenger Model, the soft supersymmetry breaking terms are de- 
cided by the gauge quantum numbers of the fields. As we will demonstrate later, this 
significantly alters the hierarchy within the neutrino mass states. In particular, we find 
that the rao dominates over m loop in case of the A' couplings but the situation is reversed 
when the R violation occurs through A couplings. This feature is characteristic of the gauge 
mediated scenario and is quiet distinct from all the earlier studies. 

We discuss the basic formalism in the next section which also contains analysis of the 
effect of the trilinear A' couplings. The third section has detailed study of the A couplings 
and we end with a discussion in the last section. 

II. FORMALISM 

We consider the following trilinear interaction in this section: 

W np = KjkLiQjdt , (5) 

where i,j, k are generation indices. In spite of very large number of these couplings, one 
could determine the neutrino masses and mixing in terms of small number of parameters if 
one assumes that all the trilinear couplings are similar in magnitude. The basic formalism 
was developed in [|HJ and we recapitulate here the relevant parts. 

The neutrinos obtain their masses from two different contributions in this case. The A' 
couplings generate radiative masses through exchange of the down squarks at the 1-loop 
level. In addition to this, the trilinear interactions also radiatively generate soft SUSY 
breaking terms which are linear in the sneutrino fields. These terms lead to additional 
contribution to neutrino masses which can dominate over the the first contribution. The 
second contribution follows from the RG improved scalar potential f[5,16 : 



V scalar = ml | t> t \ z +mi | HI Y +mj l2 | H» \* + mi H v\H\ 



,2 ~*u0 



1 



-fM B,H»H» - B^Et + H.c\ + -(gf + g 2 2 )(\ H° \ 2 - | H» \ 2 ) 2 + .... , (6) 

where we have used the standard notation for the SM fields and their masses, with B f and 
vn? vH representing the bilinear lepton number violating soft terms. Minimization of the 
above potential leads to the sneutrino vevs: 



B e .vo — m 2 tt v\ 
m 2 Li + \M 2 z cos2/? ' K ' 

where v\ and v% stand for the vevs of the Higgs fields H® and H® respectively []. 

These vevs vanish at a high scale since we are assuming only trilinear L violating inter- 
actions. They however get generated at the weak scale. The magnitudes of the parameters 
B e . and m 2 _ Hi and hence the sneutrino vevs at the weak scale are determined by solving the 
renormalization group (RG) equations satisfied by them. These RG equations are presented 
in Appendix I. The general solution of these equations can be parameterized as 



nr 



2 _ \i uD.i 



ViH x — Kpphp K ip , (8) 

where k, k' are dependent on the soft terms appearing in the RHS of the respective RG 
equations and h D are down type quark yukawa. 

The sneutrino vevs break R parity and lead to mixing of neutrinos with neutralinos. This 
in turn leads to neutrino masses. For small sneutrino vevs, the neutrino mass matrix follows 



from the seesaw approximation and is given by [20]: 

M o _ Meg 2 + g' 2 ) <v % > <vj> 

13 2(- C/i M 2 + 2M^(c + tan^)) ' [) 

where c = 5g' 2 /3g 2 ~ 0.5, M\y is the W boson mass, // is the mass term for 
Higgs/Higgsino, M2 is mass term for a gaugino, and the Weinberg angle is represented 
by 6w Assuming generation independence of the terms k, k' which was found to be a very 
good approximation in [pi, we can rewrite the above mass matrix as 

M% = m K pp h D p X' jmm h D m , (10) 



x These sneutrino vevs are derived from the tree level scalar potential. Corrections from the 



one-loop effective potential can significantly shift these naive tree level values [17]. For neutrino 
phenomenology these corrections would be important in regions in the parameter space where two 
contributions to the sneutrino vev cancel each other [fTi)] . Such regions are not encountered in 
MMM parameter space in which we are interested. Moreover, we are approximately including the 
effect of 1-loop corrections by dynamically choosing soft parameters at appropriate scale in the 
manner discussed in Refs. [17|,19[]. 



6 



where the parameters p and to are summed over the three generations and tuq now contains 
the dependence of the tree level mass on the soft SUSY breaking parameters. Only one 
neutrino attains mass through this mechanism. The other neutrinos attain mass at the 1- 
loop level. The complete 1-loop structure of the neutrino masses has been discussed in |18 



In the present case, the most dominant contributions are from diagrams having A' couplings 
at both the vertices. The mass matrix generated by these diagrams is given by 



M l tj = j^Kik^jki vi h k sin^ cos0j ln-^f . (11) 



In the above, sin 0/ cos <pi determines the mixing of the squark-antisquark pairs and M\ x 
and M^i represent the eigenvalues of the standard 2x2 mass matrix of the down squark 
system. The indices I and k are summed over. As the mixing sin0/ cos0; is proportional 
to hf , we rewrite the 1-loop contribution as, 

Mij = m loop \' ilk X' jkl h k h t , (12) 

where m loop is independent of the R violation and solely depends on the MSSM parameters. 
The total neutrino mass matrix is given by, 

M U = M° + M l (13) 

which can be rewritten in the following form when Oihf , hf ) terms are neglected: 

Mij ~ (m + m loop ) aidj + m loop hf hf A tj (14) 

where at = \' ipp hp (p summed over generations) and 

Aij = ^i23^j32 + ^32^23 — ^i22^j33 — -^J33^j22 ■ (15) 

To derive the eigenvalues of the total matrix M u , we recognise that a) The first matrix 
on the RHS of eq.(|nj) has only one non-zero eigenvalue ; b) The dominant terms in the 
total matrix M v of Oihf 2 ) are present only in the first matrix . Moreover, as we will 
show below m loop <^ mo in the Minimal Messenger Model in the purely A' case. Hence 
approximate eigenvalues can be derived up to 0( — !" s> "' oop — r), neglecting the high order 
corrections. The detailed derivation of the eigenvalues and the mixing matrix has been 



presented in |TJ|]. These eigenvalues are given as, 

m Vl ~ m loop /if hf Si 
m U2 ~ m loop hfhf 5 2 

3 

m U3 ~ (to + m loop ) J2 ' 



(m + m loop )hfJ2X' l33 2 (16) 



where 



Si = (c\ A' n - 2c 1 s 1 A' l2 + s\A' 22 ) , 
5 2 = (si A' u + 2c 1 s 1 A' 12 + c\A' 22 ) . 

The entries A',, are the elements of the matrix A' = UT,AUy where 



Uy 



( \ 

C 2 S2C3 S2S3 

-52 C2C3 C2S3 



\ -S3 C 3 y 



with 



S-2 



a 1 



a\ + a\ 



5.3 



\ 



a\ + a 2 + a§ 



The total mixing is given as ||15|| , 



K' = Uy U'y 

I 



C1C2 - S1S2C3 S1C2 + C1S2C3 S2S3 



— S 2 Ci — S1C2C3 — S\S 2 + C1C2C3 C2S3 



V S1S3 — S3C1 

where the 1-2 mixing angle 9\ is given by, 



c 3 



tan 20i 



2A' 



12 



4' — 4' 



(17) 



(18) 



(19) 



(20) 



(21) 



From eq.(|l"5D, we see that in the limit of exact degeneracy of the A' couplings, the 
parameters Aij would be zero. In this case, only one neutrino becomes massive in spite of 
the inclusion of the loop corrections. The 1-2 mixing also remains undetermined in this case. 
There is no reason a priori for the exact equality of A' and non-zero but similar value for 
these parameters determine the 1-2 mixing to be large (see eq.([H])) and also generates mass 
for the other two neutrinos. 



A. MMM and neutrino anomalies 



The parameters m loop and m appearing in eq. fll6| ) are independent of the details of the 
R violation and get determined by the soft SUSY breaking terms. We assume throughout 
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that SUSY breaking is mediated by the standard gauge interactions [gT]. We work in the 
so-called minimal messenger model P3| . It is characterized by a messenger sector with a 
pair of superfields which transform vector-like under a gauge group chosen to be SU(5) for 
unification purposes. SUSY breaking is characterized by a singlet chiral superfield whose 
scalar and the auxiliary components acquire vevs breaking supersymmetry. This breaking is 
communicated to the visible sector by loop diagrams. The gauginos acquire masses at the 
1-loop level which are given as, 

Mi(X) = &i(X) A g(x) , (22) 

where X = A < S > is the supersymmetric mass of the scalar and fermionic components 
of the singlet superfield, A is the ratio -£pc, Fs being the vev of the auxiliary field of the 
singlet. The parameter x is defined as -4. The scalars acquire masses at the two loop level. 
They are given by 

ml = 2A 2 (cial(X) + C\a\{X) + jjlf a 2 (X)) f{x) . (23) 

where 63,62 are the quadratic casimirs of the gauge groups 577(3) and 577(2) respectively 
and Yi being the hypercharge of the scalar field i , with i = {Qj, Dj, Uj, Lj, Ej, Hi, H 2 }, 
where j = 1,2,3 is the generation index. The functions f(x),g(x) are given in [23| and the 



dependence of the soft masses on these functions is minimal. In the present analysis, we 
follow [jn| and choose x — |. Since the dependence of the soft masses on x is minimal, a 
different choice of x would not significantly modify the results presented here. 

The major feature characterizing the model is the absence of A -terms and the B terms 
in the soft potential at the scale X. 

A(X) = , B(X) = . (24) 

Thus, the entire soft spectrum of this model gets essentially determined by one parameter 
A. The parameters tan/? and \i are fixed at the weak scale by requiring the breaking of the 
SU(2) x £7(1) symmetry. The relevant equations following from the tree level potential are 
given by 

2ZV 
m ik + m H 2 + 2 l J ' 2 

r, m 2 fr — m 2 TT tan 2 8 1 . . 

^ = -Jh — -^ tL _ _M\ (25) 



sin 20 



where all the parameters on the RHS of the above equations are evolved to the weak scale 
using the MSSM RGE. Because of the boundary condition eq.(^), the value of the B 

9 



parameter at the weak scale remains small. This pushes tan/3 to very large values in this 
model ||T§U25| . The above equations are strictly valid in case of the MSSM and R violating 
trilinear couplings can give corrections to these. The smallness of neutrino masses however 
require very tiny trilinear couplings of O(10~ 4 — 1CT 5 ) which would contribute negligible to 
eqs.( p5f) . We thus continue to use eqs.([25|). 

The neutrino masses in eq.(|i|) are strongly hierarchical in the limit m loop <C m . Specif- 
ically, one obtains from eq.fll6|), 



Wi nnn Tfl 



loop "^ I ^ ] (26) 

//^3 m m b VEi(A- 33 ) 2 / 

where m s and mb represent the strange and the bottom quark masses. For all the A' of 
similar magnitudes, this ratio is completely determined by the parameter A. 



We can determine the above ratio by exactly solving the RG equations, (p5|) . Before 
doing this, it is instructive to study the approximate expressions obtained when one neglects 
the Q 2 dependence of the parameters appearing on the RHS of the RGE. In other words, 
we neglect the effect of running of the soft masses (from high scale, X to weak scale M£) 
appearing in the expressions of the RGE as well as for neutrino masses. Instead, we take 
them to be their high scale values given by eqs.(^). Noting that more dominant contribution 
to sneutrino vev comes from the m 2 _ Hi term in eq.(0) and integrating the corresponding RGE 
for m 2 ,. Hl in the above approximation, one finds in this simplifying case : 



3tt 2 J A a 5 2 (X) \ M 



The m loop defined in eq. (|TTJ) has the following approximate form in the same approximation 
as above where we neglect the running of the soft masses. 

(v 2 a\ cos/? sin/? 
mi - ~ { -W) 8 n 2 al(X) (28) 



Eq . (|27|) clearly demonstrates that the often neglected [|13| , |T4! RG induced contribution 



dominates over the loop contribution in the present case of MMM just as in the case of 



the supergravity induced breaking pT5| , p6| . One would have naively thought that this will 
not be the case in gauge mediated model since the running of masses in this case (signified 
by t = ln(-^-) ~ 11 in eq.(p7|)) is over much smaller range than in the supergravity case 
where t = ln( -fjM 21 ) ~ 66. But smallness of t in MMM is compensated by the largeness of 
the ratio —%■ (signified by the factor) -$ . As a result of which the value of m here can be 



comparable to the corresponding value |L5| in supergravity case. From the expressions we 
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see that dependence on the A is more severe for the 1-loop mass, m loop compared to the tree 
level contribution, m . However, the /i parameter in the numerator increases approximately 
linearly with A [19"]. This makes the A dependence of the both the contributions essentially 



the same. The ratio m ' oop following from eqs.(^,^) is given as: 

m loop ( 7t2 \ ( A 4 tan/5 v 2 \ ( a%(X)\ 



(29) 



m V 3 / V M w A t 2 J V«3P0/ 

From the above we see that the dependence of the ratio on A is essentially determined by 
the way the \x parameter scales with respect to A. This leads to a very mild dependence of 
the ratio on A. Such A 'independence' has also been seen in the case of bilinear R-violating 
models in MMM [§. For A = 100 TeV, t = 2 In (§-") = 10.6, tan/3 = 46 and /i = 400 
GeV, from the above we see that the ratio is 0.39. From eqs.( |2~7j , |2"8| ) we see that the typical 



order of magnitude for the ratio of the mass eigenvalues, eq.(|2"6|) is : 

■ 10~ 3 (30) 



ITlm , — 3 



m U3 



However the above expressions are approximate. We have determined this ratio exactly by 
solving the relevant RG equations numerically. For this, we closely follow the work of [H5| 



where two loop RGE for the B parameter were used for fixing the sign of fi parameter at the 
weak scale. As mentioned above, we have chosen x = ~. This fixes X = 2A. Following [19 



the decoupling scale 'Q ' is defined to be the geometrical mean of mg(A),m^(A). Such a 
choice of Qo would be approximately equivalent to adding corrections due to the complete 
1-loop effective potential. The gauge couplings and the Yukawa couplings are evolved from 
M z to Qo using the Standard Model (SM) beta-functions and using MSSM beta-functions 
from Qo to the high scale, X. All the soft masses and soft parameters defined at X are 
evolved from X to the decoupling scale Qq. The ratio — ^, determined following the above 
procedure is plotted in Figure 1 for A varying from (50 — 150) TeV. 
From the figure we see that the ratio of the eigenvalues, 



is typically around the expected value, eg. (|30|) . While the ratio shown in Fig. 1 is completely 
fixed by the value of A, the neutrino mass ratio is uncertain by a number of 0(1) which is 
related to the trilinear parameters. 

We now turn to discussing feasibility of the model for the simultaneous description of 
the solar and atmospheric data. The two generation analysis of each of these experiments 
constrain the value of the relevant (mass) 2 difference and mixing angle. The allowed values 



for the parameters p8H30|] are displayed below: 
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Anomaly 


Solution 


Am 2 {eV 2 ) 


tan 2 6» 


Solar 


MSW-SMA 


(2- 10) x 10~ 6 


(1-20) x 10" 4 




MSW-LMA 


(2 - 80) x 10~ 5 


0.2 - 4. 




LOW-QVO 


1 x 10- 10 -4 x 10~ 7 


0.1-8. 




Vacuum (Just- So) 


(4 - 12) x 10~ 12 


0.1 -7. 


Anomaly 




Am 2 (eV 2 ) 


sin 2 29 


Atmosphere 




(1-8) x 10~ 3 


0.83- 1. 



At present, global analysis including the recent results from the day /night recoil electron 

energy spectrum and charged current rates from SNO of the solar neutrinos from super- 

Kamiokande [[| favours the MSW-LMA solution but the other solutions are also allowed 

28|,30]. It follows from eq.fl3~l|) and the table that the most natural solution for the solar 



neutrino problem is through the vacuum oscillations but quasi-vacuum solution can also be 
obtained, if the A' dependent factor in eq. (j3lD is somewhat large (e.g. ~ 5) instead of exactly 
being one. The most preferred LMA solution cannot however be obtained. 

The mixing among neutrinos is essentially controlled by ratios of trilinear couplings. 
This mixing is given by eq. (|20"D . It is seen from this equation that a choice of angles .§1,2,3 is 
possible which reproduces two large and one small mixing angles as required by the present 
data. As an example, consider the choice 



c\ = c 3 = si = s 3 =-7= ;s 2 ~.13 



(32) 



This gives 



w?20 A = *K'J{\-K 



L /x3 



0.99 



sin 2 26 s = AK e 2 K e2 2 ^ 0.95 



Kb = 0-09 • 



(33) 



This choice reproduces the required mixing angles and also satisfies the CHOOZ constraint. 
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III. MODELS WITH \ IJK 

In this section, we discuss the structure of neutrino masses and mixing in the presence of 
only trilinear A couplings. The lepton number violating part of the superpotential is given 
as, 

W$ = K^Lpl . (34) 

There are two basic changes here compared to the last section. Firstly, the A^fe are antisym- 
metric in the first two indices restricting their total number to nine. This strongly restricts 
neutrino mass structure and one does not get phenomenologically consistent spectrum when 
all the trilinear couplings are assumed to be similar in magnitudes. Secondly, unlike in the 
A' case, the loop induced contribution dominates over the tree level SUSY breaking in the 
minimal messenger model. Such dominant loop contribution has been earlier seen in some 



particular regions of the parameter space in the mSUGRA framework [|15|,0]. Here this 
dominance follows generically. 

We present here a formalism which is similar to the spirit of the previous section in order 
to understand the basic features of the neutrino mass matrix . As before, the RG improved 
effective soft potential contains terms B ti and rr? v . Hx which break the lepton number. These 
terms are generated due to the presence of A couplings in the superpotential. At the weak 
scale, the magnitude of these depends on the respective RGE, which we have presented in 
Appendix I. The solutions of eqs. (|66"D can be written as, 

-D ei = AipptlppKip, l"J^J 

m u l H 1 = ^ipphppKip , (36) 

where (i ^ j) due to the anti-symmetric nature of the A couplings and k and k' represent the 
dependence on the soft masses in the RGE and h E are the charged lepton yukawa. Following 
similar arguments as in Section II, the presence of these terms in the scalar potential would 
lead to a tree level neutrino mass matrix of the following form: 

M% = mAh , (37) 

where mo contains the dependence on the soft terms and hi are given as, hi = \i VV h E (i ^ p). 
In addition to the tree level mass, the presence of A couplings also gives rise to contributions 
at the 1-loop level. Assuming only canonical 1-loop contributions to be the most dominant 
contributions the 1-loop level mass matrix has the form: 
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1 M 2 

•M;„- = ^-^X Uk X jk i v 1 h k sulfa cos (pi In 



M,V 



(38) 



lJ 167T 2 

where sin 0; cos 0; and Mu,M 2 i represent the mixing and the eigenvalues respectively of 



the standard 2x2 stau slepton mass matrix. Following the previous section, we rewrite the 
above as, 



■M-ij — Tiioop XukXjkihi h k 



(39) 



In writing the above, we have implicitly assumed the anti-symmetric nature of the cou- 
plings. The total neutrino mass matrix is given as, 



M = M [) + M 



i KAl 



The above can be rewritten in the form: 



,£ 2 uEuE\ 



Mij = (m + m loop )bibj + m loop hfh^B i:j + 0(h% , hfhf), 



(40) 



(41) 



,£ 2 uEuE\ 



where we have neglected ©(/if , h 2 hf) contributions to the mass matrix. The matrix B is 
given as 



A132A123 — A133A122 A123A232 — A122A233 A132A323 — A133A322 



B 



A103A 



123^232 — ^233^122 



A 933 A 1 



A132A323 



A 133 A3 



















(42) 



\ ^132^323 — ^133^322 " " / 

We diagonalise the total matrix Ai in the same manner as for the A' case and as described 
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15| . However we do not make any assumption on the relative magnitude of mo and m loop . 
The approximate eigenvalues correct up to 0{h 2 2 ) are derived as, 



m ui "-- 



m loo XKvt + OQ4 



E 2\ 



m U2 « 



E 2\ 



m loop h%h2r) 2 + 0(h h 2 
3 
m U3 as (m + m loop ) ^ b 2 

i 

« (m + m loop )/if % + 0{h$h%) 



(43) 



The parameters rji,r]2,rj 3 are given by, 



m 



{c\ B' u - 2c lSl B[ 2 + s\B' 2 



22 j 



(s 2 1 B' n + 2c 1 s 1 B[ 2 + c 2 1 B' 22 ) 



2 

?33 



(44) 



j=l,2 
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The parameters B\, are elements of the matrix, B' = U^BU\ where, U\ has the same 



form as U\> in eq.(|T^) with the angles now given by, 
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Si 



bj + bj 



S3 



b\ + bl 



\b\ + bl + bl 



(45) 



The total mixing matrix is given as in eq.(^0|) with the 1 — 2 mixing angle given by, 



tan 20i 



2B' 



12 



-B22 - B' u 



(46) 



In the above, we have made no assumption on the relative magnitudes of mo and m loop . 
Only assumption is almost equality of all the A couplings. Most definitions are formally 
the same as in case of the A' couplings but the physical consequences are quite different. 
This follows from the expressions of mixing angles ^2,3 in eq.(^). If no hierarchy is assumed 
among the A couplings then due to antisymmetry of A we get, 
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h 



0( 



hf 2 ' 



As a result, eq.(H3J) implies C3 ~ 0(^-) and the mixing matrix, eq.(|2"0|) can be written 



as 



K = U X U' X 
I 



\ 



\ 



,m, 



m T 



(47) 



/ 



CiC 2 SiC 2 s 2 
-S 2 Ci — SiS 2 c 2 

si c\ 
The equations for 771,772 and tan20i take the following approximate forms in this case: 

771 fa s\clBu + 2.6120251(01 - sis 2 ) 
7] 2 ^ c\clBu - 2Bi 2 CiS 2 (s 1 + CiS 2 



(48) 



tan 20i 



-2B- 



12 



(49) 



B u c 2 - 2s 2 B l2 

With the hierarchical masses, the effective mixing angles 9a and Oqhooz probed by the 
atmospheric data and the CHOOZ experiment respectively are given by 



sin 2 29 A = AK%{1 - K%) sin 2 29 C hooz « 4K 2 3 (1 - K 2 e3 ) . 

Eq . (|47|) implies that these two mixing angles are equal in conflict with the observation. 
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The above derivation has not assumed any specific mechanism for the supersymmetry 
breaking and thus the conclusions are valid in supergravity scenario as well as in the MMM 
case considered here. It is quite interesting that in spite of the presence of nine indepen- 
dent A parameters, one cannot explain solar and atmospheric neutrino anomaly as long as 
these parameters are of similar magnitudes. Departure from equalities of A can lead to 
explanations of these anomalies and we shall present a specific example in the next section. 

A. MMM and Neutrino Anomalies 

In this sub-section, we first determine the numerical values of the parameters m loop and 
mo which enter the neutrino mass. The dependence on trilinear couplings is factored out in 
defining these parameters. But their numerical values are quite different here compared to 
the A' case studied in the earlier section. This follows since the soft sfermion masses which 
determine these parameters are quite different in these two cases ( see RG equations in the 
appendix ). 

In models with gauge mediated SUSY breaking, the soft masses are proportional to the 
gauge quantum numbers they carry. Thus particles with strong interactions have much larger 



soft masses compared to the weakly interacting particles, as is evident from the eqs.(22,P3|). 
The effect of the gauge couplings also trickles down to the parameters B e .,m 2 ,. Hi through 
the corresponding soft masses present in their respective renormalization group equations. 
In the presence of purely A' interactions, the strong coupling determines the magnitudes of 
B ti , rnl iHl at the weak scale and in turn the tree level mass as is evident from eg. (|27|) . The 
loop contribution is still however determined by the weak coupling, eg. (|28|) . It is thus the 
interplay between the strong coupling and the weak coupling which lead to a suppressed 
loop mass in the case of purely A' couplings. In the case of pure A couplings, the squark 
masses do not enter the definition of mo and m loop and both of these are determined by 
the weak coupling. However the dependence on the power of the weak coupling is different. 
Following the same method as described in the A' case (above eq.(pTD), an estimate of the 
parameter itlq is given by, 

(cos0\ 2 M? v 1 / X'V 

m »=U^j — HDrMl) (50) 

The 1-loop contribution m loop can also be estimated in the similar manner as, 



_ fv 2 fi\ cos/? sin g 
m loop - ^— 24 ^ (x) (51) 
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The ratio m ' oop is then given by, 



"'o 



m loop /8?r 2 \ / V // tan/3 \ 1 , , 

m ~ ^ 3 J ^ £ 2 A M^ J a 2 (X) [ } 

where t = In Qg-Y The above is typically of O(10 2 ) for A = 100 TeV, /i = 400 GeV, 
tan /3 — 46 which shows that the tree level mass is much suppressed compared to the 1-loop 
mass. Comparing the above equation with that of the corresponding one for the A' case, 
eq. (|29"D , we see that the absence of strong- weak interplay in this case leads to a much larger 
ratio. 

In a general mSUGRA inspired scenario, the tree level mass is much larger compared 



to the 1-loop mass |15[ for large range in MSSM parameters, irrespective of the nature of 
R parity breaking. In the present case, the relative importance of loop and the tree level 
contributions is sensitive to the nature of R violation as demonstrated above. This feature 
arises not as a consequence of running of soft masses but due to difference in the relevant 
RG equations in case of A and A' couplings and the boundary conditions themselves which 
strongly depend on the gauge couplings in these models. 

We have determined the ratio m ' oop by solving the relevant RG equations numerically in 
the manner described in section (3). This ratio is plotted versus A in Fig.2 for A = (50 — 150) 
TeV. In this range, 

m ioop 



25-45. (53) 

m 

as expected from eq. (|52"D . 

The dominance of m loop has the following important implication. The neutrino mass 
ratio following from eq. (|43[) is given by 

m U2 _ m loop rrif, rj 2 , , 

m Vz m + m loop m T r] 3 

In the previous case, the hierarchy in m loop and mo resulted in strong hierarchy between 
neutrino masses. As a result, one could only obtain vacuum or quasi- vacuum solution for 
the solar neutrino. Here due to m loop ^> mo, hierarchy in neutrino mass is much weaker, 
and we have, 



j i 



^^~5x 10- 2 (55) 

m v , m T 



We can thus easily get the scale relevant for the LMA solution of the solar neutrino. As 
already argued the mixing pattern is not appropriate if all the A couplings are similar. This 
is no longer true if A obey some specific hierarchy as we discuss now. 
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B. Illustrative Model 



We assume that the couplings A123, A233, A322 dominate over the rest and neglect the 
latter. Moreover we assume that non-zero couplings satisfy the following hierarchy, 

A 233 „fh£\ *-- /ur 



t<-°& 



(56) 



A322 \h 3 J A322 

We do not have strong theoretical reasons to assume the above hierarchy The follow- 
ing considerations should therefore be viewed as an example which leads to the successful 
explanation of the neutrino anomalies. 

The tree level mass matrix in the presence of these couplings is given by, 



/ 



M[ 



m 












\ 2 h E 2 
A 233' i 3 







\ 



A233^-3 -^322^2 



Ao.SS/iq A322^9 Aooa/ln 



^233^3 ^322"-2 

The 1-loop level mass matrix is given by 



v 322"-2 



(57) 



Ml 



711 



loop 







Al23^2 A232^3 



Al23^2 -^232^3 -^233^3 A232^2 -^323^3 



V 



A 2 32^2 -^323^3 ^322^2 



(58) 



/ 



In view of the hierarchy in eq. (p5]) , the total mass matrix has the following simple form 



M* 



( 




X 


\ 




X 


A 


A 


y° 


A 


A » 



(59) 



where 



A = (m + m loop )\ 2 233 hf 2 , 
x = —m loop Xs22^i23h 2 h 3 . 



(60) 



We shall assume x < A which is consistent with the hierarchy in eq. ([56]) . One can diagonalise 
the above matrix: 



R 



R^ 



12^12j ^131^13 



#23(tt/4) M s [R 12 (6 l2 ) R 13 (6 13 ) i? 23 (7r/4)] J « Diag.(m Ul ,m U2 ,m U3 ) 



Here, Rij denotes rotation in the ij plane with angle %. We have neglected a small 
contribution of 0( % ) to the 2 — 3 mixing angle in the above derivation. The mixing 
angles are given by 

tan26> 12 ~ , 

x 

tan26n w ~ij2 ■ (61) 

The eigenvalues can be approximately written as 

o 
X X 



m U2 w 



2 

x ar 

72~8A 
;r 2 



m„ 3 n2A+—. (62) 

It is seen from the last two equations that all the mixing angles and the masses are 
predicted in terms of only two parameters namely, x and A. The atmospheric mixing angle 
is predicted to be around ir/4 and the other two mixing angles can be expressed in terms of 
the solar and atmospheric scales. Using eqs. (|6l|j62]) we find, 



f oa 4y/2A /- /Sv^A^A 1/3 
tan 29 80 i ar w — — w 4V2 I I 

tan2^ 00 ^ I - 7 =.- 7 =^^^j . (63) 

Choosing A A ~ 3 • 10- 3 eV 2 and A soiar ~ (2 - 50) • 10- 6 eV 2 we get, 

j* (0.19 -0.57), 

tan 2 6 solar w (0.93 - 0.81) , 
U e3 (CHOOZ) w (0.06 - 0.19) . (64) 

The predictions for the mixing angles are in very good agreement with the observations 
which prefer large mixing angle solution for the solar neutrino. The required value of ^ is 
also consistent with the assumed hierarchy in eq.(|56|) among the trilinear couplings. 



IV. DISCUSSION 

We have discussed the structure of neutrino masses and mixing in the Minimal mes- 
senger model (MMM) of gauge mediated supersymmetric breaking with purely trilinear R 
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violating interactions. We considered two specific cases of purely A' interactions and purely 
A interactions for simplicity. The model contains very large number of parameters even 
under this simplifying assumptions. Remarkably, it is possible to make meaningful state- 
ment on the neutrino spectrum in spite of the presence of many unknown parameters if all 
these parameters are assumed similar in magnitude. This is a natural assumption in the 
absence of any specific symmetry to restrict the trilinear R parity violation. It is not always 
easy to justify this specific choice, e.g use of a U(l) symmetry which uses Froggatt -Nielsen 
[FN] mechanism to obtain quark and lepton masses tend to forbid all the trilinear terms 
altogether |Q. 

In the case where only A' couplings are present, one naturally gets large mixing between 
the neutrino states. Further, the MMM offers a very constrained structure giving rise to 
a large hierarchy between the masses ~ O(10~ 2 ) for all the parameter space. The model 
is suitable for obtaining simultaneously solutions for atmospheric neutrino problem and 
quasi-vacuum oscillations. 

Assumption of approximate equality of A couplings in case with only A couplings, leads 
to very constrained and phenomenologically inconsistent pattern for neutrino mixing. This 
conclusion follows on general grounds and it is true even if SUSY breaking is induced by 
supergravity interactions. It is quite interesting that one can arrive at this strong conclusions 
in spite of the presence of many unknown parameters by simply assuming them to be of 
similar magnitude. 

One can obtain consistent picture of neutrino anomalies if A couplings are assumed to 
be hierarchical. We provided an example which leads to two large and one small mixing and 
correct hierarchy between the solar and atmospheric neutrino scales. 

One interesting result of this analysis is the interplay between the sneutrino vev induced 
contribution and the loop induced contribution to neutrino masses. In the context of su- 
pergravity induced SUSY breaking, it has been shown that the large logarithmic factors 
induced due to RG scaling enhance the sneutrino vev induced contribution compared to the 
loop contribution. We showed that this remains true even in the gauge mediated models of 
SUSY breaking in case of the trilinear A' couplings. In the mSUGRA model, the dominance 
of tree level mass follows simply from the large factor t = In ffi T in sneutrino vev gener- 
ated by running of soft parameters. In the present case, the tree level dominance occurs 
essentially due to boundary conditions. In case of A' couplings, m is determined by squark 
masses which depend upon a^{X). mo dominates over loop contribution in this case. For 
A couplings, mo is determined by sleptons rather than by squarks masses. Due to their 
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dependence on weak couplings, slepton masses are much smaller than squark masses. As a 
consequence, m is suppressed compared to A' case. This results in loop dominance if R is 
violated by A couplings. 

The neutrino mass hierarchy strongly depends on the ratio m ' oop . In case of the tree 
level dominance (purely A' couplings) one obtains strong hierarchy and vacuum solution 
while the case with loop mass dominating corresponds to milder hierarchy and the LMA 
MSW solution. This analysis along with other similar analysis |T5|J16|j2^| therefore underlines 
the need of including both the contributions to neutrino masses in a proper way. 
Acknowledgments: We warmly acknowledge discussions with U. Chattopadhyay which 
have been very helpful. 

Appendix I 

Here we present the Renormalization Group Equations (RGE) for the soft parameters B e . 
and m 2 . H for the two cases considered in this work: either purely A' couplings or purely A 
couplings are the sources of lepton number violation in the superpotential. These equations 
have been derived using the general formulae given in [[Jl[] . These equations can also be found 
in [|16[], whereas the equations for the standard soft parameters appearing in the RHS of the 
equations can be found in many papers such as [EE3]. In writing the below eqs. flB5|) (eqs. flSED ) , 
we have neglected 0(X 2 ) (0(X 2 )) corrections. The notation is as described in the text. 
A' couplings in the superpotential: 

dBeM B ei (t) (~Y t (t) - hf® + h 2 {t) + iai(«) 



dt lV ' V 2 w 2 l w 2 w 10 

' Kt)KiMhUt)(lBJt)+A}'(t) 



d ^f^ L = < Hl (t) {-\y, E it) - \Y b (t) - \Y T (t)) - ^- 2 K PP (t)h d pp (t) «(t) 

+ m 2 Li (t) + 2 m 2 Qp (t) + 2 A^{t)A^{t) + 2 m 2 Dp (t)) (65) 

A couplings in the superpotential: 

dBeXt) ~ BeM (-hit) - l -Yf(t) + h 2 (t) + iai(t) 



dt lV ' V 2 w 2 l w 2 w 2 

V(t)\ ldd (t)h% d (t) (A^t) + l-B^t) 



Jfcr2 rv '"™ ,v "'"«« ( "' \--iaa\-; ■ c^ 

^f^ = < Hl {t) (-\Y % E {t) - \Y T (t) - \Y b {t)) - ^X^mm (m 2 Hl (t) 

+ m 2 Li {t) + 2 m%{t) + 2 A^{t)A%{t) + 2 m%(t)) (66) 
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FIG. 1. Neutrino mass eigen value ratio —^- plotted versus A (GeV) assuming trilinear A' 
couplings of similar strength to be the only source of R violation. A is defined in the text. 
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Fig. 2 
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FIG. 2. The neutrino mass ratio 



m-ioo 



plotted versus A (GeV) assuming trilinear A couplings of 



similar strength to be the only source of R violation. 
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